COUNTING NUMERICAL SETS 
WITH NO SMALL ATOMS 



JEREMY MARZUOLA AND ANDY MILLER 

Abstract. A numerical set S with Frobenius number g is a set of integers 
with min(5) = and max(Z — S) = g, and its atom monoid is A{S) = 
{nSZ|n + sSS for all s S S}. Let ■yg be the number of numerical sets S 
having A{S) = {0} U (g, oo) divided by the total number of numerical sets 
with Frobenius number g. We show that the sequence {79} is decreasing and 
converges to a number -yoo ~ .4844 (with accuracy to within .0050). We also 
examine the singularities of the generating function for {"fg}. Parallel results 
are obtained for the ratio 'yg of the number of symmetric numerical sets S with 
A(S) = {0}U(g, 00) by the number of symmetric numerical sets with Frobenius 
number g. These results yield information regarding the asymptotic behavior 
of the number of finite additive 2-bases. 



Let Z denote the additive group of integers and let N denote the monoid of 
nonnegative integers. Both of these sets are hnearly ordered by the Archimedean 
ordering and we will use standard interval notation to describe their convex subsets. 
If n S Z and S C Z then the translate of S* by n is the set n + S^{n + s\ sE S}. 

A numerical set S" is a cofinite subset of N which contains 0, and its Frobenius 
number is the maximal element in the complement N— 5* Equivalently, a numerical 
set S with Frobenius number g is a set of integers with min(5) = and max(Z— 5*) = 
g. A numerical set which is closed under addition is called a numerical monoid. 
Every numerical set S has an associated atom monoid A(S) defined by 

A{S)^{neZ\n+SCS}, 

and this is easily seen to be a numerical monoid with the same Frobenius number as 
S. Note that A{S) C S and that S" is a numerical monoid if and only if A{S) = S. 
The nonzero elements of A(S) are referred to as the atoms of S. 

For each g > let Ng be the numerical monoid 

N<,-N-[l,.g] = {0}U(g,oo) , 

which has Frobenius number g when g > oj^ The atom monoid of every numerical 
set S with Frobenius number g contains Ng and the complement 5 — Ng is a subset 
of (0,g). Conversely, the union of Ng with any subset of {0,g) is a numerical set 
with Frobenius number g. Therefore the set 

S{g) = {S" C N I 5 is a numerical set with Frobenius number g} 



^This definition differs from that employed in |AM| where a 'numerical set' would be a translate 
n + S of a numerical set S (in the sense given here) by an arbitrary integer n. Since the atom 
monoid of n + S equals the atom monoid of S, this variation of the definition should not lead to 
any confusion. 

■^The Frobenius number of No = N is — 1, and this is the only numerical set with nonpositive 
Frobenius number. 
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is in one-to-one correspondence with the power set 7^(0, g) consisting of all subsets 
of (0, g), and S{g) has cardinality 2^^^. The subset of S{g) consisting of numerical 
monoids is a much more difficult set to enumerate. This is examined in Backelin's 
paper [B] where it is shown that for large values of g roug hly 3 X 2L(ff-i)/2J of the 
2^~^ elements of S{g) are numerical monoids. 

If M e S{g) is a numerical monoid then the anti-atom set of M is the set 

g{M) = {Se S{g) I AiS) = M} . 

This is contained in the larger set S{M) = {5 G S{g) \ M C ^4(5*)} whose elements 
might be considered to be 'Af- modules']^ Notice that S{g) — S{Ng) and we will 
also write g{g) = g{Ng). This paper is motivated by the following question which 
we shall refer to as the Anti-Atom Problem. 

For a given numerical monoid M with Frobenius number g 
how many numerical sets in S{g) have atom monoid M? 

Thus, for a given monoid M , the Anti-Atom Problem asks to compute the cardi- 
nality of Q(M). This problem is certainly unwieldy given that it fundamentally 
presupposes an enumeration of the set of numerical monoids in S{g) — an enumer- 
ation which Backelin has shown to be intractable at best. Nevertheless we will be 
able to frame aspects of the problem in a clearer light. Our main result will show 
that there is one monoid M in S{g) (that monoid being AI = Ng) which itself is the 
atom monoid for approximately 48.4% of all numerical sets in S{g) for large values 
of g. In order to describe this in more depth we first need to discuss symmetry 
and pseudosymmetry in numerical sets. These concepts are important throughout 
much of the theory of numerical monoids and numerical sets (see |FGHj , |AMj and 
[X] for example), and will play a role in many of our discussions. 

A numerical set S G S{g) is symmetric if an integer x is an element of S if 
and only if g — x is not an element of S*. In other words, S is symmetric when the 
reflection on Z given by x g — x carries S onto its complement Z — 5. Notice 
that only numerical semigroups with odd Frobenius number can be symmetric. 
A numerical set with even Frobenius number g is said to be pseudosymmetric if 
g/2 ^ S and for each integer x ^ g/2, x is an element of S if and only if 5 — a: is 
not an element of S. Symmetry and pseudosymmetry can also be described using 
the notion of duality of numerical sets. If S" € S{g) then the dual of S is the 
numerical set S* ^ {n E Z \ g — n ^ S} , and it is not hard to show that S* £ S{g) 
and that A{S*) = A{S) (more background can be found in section 1 of |AM] ). The 
numerical set S is symmetric if and only if S* = S, and it is pseudosymmetric if 
and only if g is even and S* = S (J {g/2}^ For each numerical set S G S{g) there 
is a rational number type{S) no smaller than one, called the 'type of S", which 
satisfies the property that S is symmetric if and only if type{S) = 1. The type of a 
numerical monoid AI £ S{g) is always an integer, and it equals the cardinality of its 
omitted atom set 0{M) = {n £1 - M \ n + {M - {0}) CM}. Since 0{M) C N 
and g £ 0{A4), the type of a numerical monoid M £ S{g) is in the interval [i,g], 
and the largest possible value type{Ad) — g is only achieved when M = Ng. The 



''in IBF| the elements of cS(Af) are called 'relative ideals over M'. 

^ More generally, if the symmetric difference of S and S* contains no more than one ele- 
ment then S is symmetric, pseudosymmetric or "dually pseudosymmetric" (meaning that S* is 
pseudosymmetric) . 
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following elementary results allow us to solve the Anti-Atom Problem for symmetric 
and pseudosymmetric numerical monoids. 

Proposition 1. Suppose that M is a numerical monoid and that S is a numerical 
set with A{S) = M. Then M CS C M* . 

Proof. Let S* be a numerical set in S{g) with A{S) — M and s G S*. If g — s were 
an element of AI then g ~ s + {g ~ s) would be an element of S, which contradicts 
g being the Frobenius number of S. Thus g ~ s <^ M which implies that s e M* , 
and M = A{S) C S C M*. □ 

Corollary 2. A numerical monoid M €E S{g) is symmetric if and only if there is 
just one numerical set (which must be M itself) whose atom monoid is M. If M is 
a pseudosymmetric numerical monoid then there are precisely two numerical sets 
(which must be M and M*) whose atom monoid is M. 

Proof. Let M E S{g) be a monoid. If M is not symmetric then AI ^ M* but 
A{M*) = A{M) = M, and so there are at least two distinct numerical sets in Q{M). 
On the other hand, if M is symmetric and S £ Q{M) then AI C S C M* — M and 
S = Ad. If AI is pseudosymmetric and S G Q{M) then Ad <Z S <Z M* = A/I \J {g/2}, 
so that S equals M or AI* . □ 

This corollary then provides the first positive answers to the Anti-Atom Prob- 
lem: namely, that |CJ(M)| = 1 when Ad is symmetric and that — 2 when 
Ad is pseudosymmetric At the other end of the spectrum, we shall show that the 
anti-atom set of (which is the numerical monoid in S(g) farthest removed from 
being symmetric since type{Hg) = g is the largest possible type among all monoids 
in S(g)) is an order of magnitude larger in size than that of any other numerical 
monoid with Frobenius number g. To establish this we will examine the sequence 
7g = We introduce a combinatorially defined sequence of positive in- 
tegers {Afe} with the property that 1 — 7^ is a partial sum of the convergent infinite 
series X^fcLi ^fc4~*^. This allows us to show that {73} is a decreasing convergent 
sequence and that the limit 700 is approximately equal to .484451, give or take 
.0050. The integers Ak are combinatorially related to integers A'^. which turn out 



to equal (theorem 10 1 and there is a nice recursive relation between these 

two sequences (theorem 111. This relationship will enable us to obtain informa- 
tion about the singularities of the generating functions for the sequences {Ak} and 

In addition to forming a large subset of S{g), the numerical sets in Q{g) have 
nice properties in terms of the direct sum decompositions discussed in |AMj . Given 
numerical sets S and T and relatively prime atoms a e A{S) and h £ A{T) the 
direct sum of S and T is the numerical set bS © aT = {6.s + at | s G 5 and t G T}. 
Every numerical set S can be trivially described as S ~ \S ® oN for any nonzero 
a S A{S), but if this is the only kind of direct sum decomposition of S then we say 
that S is irreducible. Every numerical set can be expressed as a finite direct sum of 
irreducibles. By [AMI Proposition 4.4], the only numerical set in [J{G{g) | > 1} 
which is not irreducible is Ni = 2N © 3N. Thus our results show that at least 
47.94% of all numerical sets in S{g) are irreducible. Another nice property is that 



^We showed that |CJ{A/)| = 1 if and only if M is symmetric, but it is not hard to construct 
numerical monoids M with I^CM)! = 2 that are not pseudosymmetric. 
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the type function is multiplicative when restricted to I ^ 1} by |AMI 

Proposition 5.3] (that is, the type of a direct sum is the product of the types of its 
factors, if the factors have no small atoms). Multiplicativity of type was a central 
theme in [AM ]. We also mention that when a numerical set S is in G{g) its type 
can be computed via the formula 

_ |5n[o,g)||^-n[o,g)| 
type[b)- ^sns*n[0,g)\^ ' 

which is readily derived from, but much simpler than, the general formula for the 
type of an arbitrary numerical set given in |AMj . 

In the last two sections of the paper we explore parallel ideas for counting the 
number of symmetric numerical sets in G{g). This study is suggested by Backelin's 
examination of the number of symmetric numerical monoids in S{g) in pj. We 
show that the ratio "fg of the number of symmetric numerical sets in Q (g) by the 
total number of symmetric numerical sets with Frobenius number g has a limit 
7fJ, which is approximately equal to .23644. We also obtain information about the 
singularities of the generating function for the sequence {7^'^}. In many ways the 
analysis of {7^'^} turns out to be more elementary than that of {jg}- For example, 
in the symmetric setting we obtain two recursively related sequences of integers 
{A^} and {A^'}, and the odd terms in the second of these sequences coincides 
with a well-known sequence consisting of the numbers of additive 2-bases for k. 



Numerical sets with no small atoms 



Let S* be a numerical set with Frobenius number g. A small atom for S" is a 
(nonzero) atom for S which is less than g. A small atom for S is said to be large 
if it is greater than g/2. The first result says that every numerical set which has a 
small atom will have a large small atom. 

Lemma 3. Let S be a numerical set in S{g). If S has a small atom then S has a 
small atom larger than g/2. 

Proof. If g is even then g/2 is not an atom of S since g/2 + g/2 ^ S. Suppose S 
has an atom less than g/2 and let k be the largest such atom. Then 2k is a small 
atom of S, and 2k is greater than g/2 by the choice of k. □ 

The set S{g) is partitioned into two subsets 

Q{g) = ^(Ng) = {5 G S{g) \ S has no small atoms} 

and 

^{g) — {S E S{g) I S has at least one small atom} . 
For each g > 0, Ng G G{g) and G{g) is nonempty. On the other hand, B{g) contains 
all of the numerical monoids in S{g) other than Ng, and B{g) is nonempty when 
g > 2. We are interested in the two ratios 

^ |g(g)| ^ |g(.g)| 
\S{g)\ 29-1 

and 

^ ^ \G{g)\ ^ |g(ff)| 
\S{g)\ 29-^ ■ 

Observe that < /3g,7g < 1 and that f3g + jg = I. 
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For each S E S{2n ~ 1) and e G Z2 = {0, 1} we define 

S', = (STl [0,n- 1]) U{en}U (l + S'n [n,oo)). 



* * * 



n-1 n 



2n-1 



M+l 



n-1 n n+1 



M+l 



2n 



n-1 n n+1 



2n 



Figure 1 . The numerical sets S'q and S[ 



Lemma 4. The correspondence {S, e) i— > S'^ is a bijection from S{2n — 1) x Z2 to 
S{2n) which carries Q{2n — 1) x Z2 onto Q{2n). Furthermore, "f2n = 72n-i and 

P2n = /32ri-l- 

Proof. The correspondence (S*, e) 1— > 5*^ is injective by definition, and it is also 
surjective: if S' e S{2n) then S' = S'^ where S is the union of S' n [0,n — 1] and 
-1 + (5 n [n + 1, 00)), and e equals if n ^ 5' and 1 if n e S". 

It is not difficult to see that an integer a; is a large small atom for 5* if and only 
if 1 + a; is a large small atom for S'^. By lemma [s] this implies that a numerical set 
S € S{2n — 1) is in Q{2n — 1) if and only if 5*^ is in Q{2n). To complete the proof, 
we note that 72™ = |^;(2n)|/22"-i = 2\g{2n - l)|/22«-i -i2n-i- □ 

For integers g and k with g > fc > 0, let 

B{g, k) = {S E S{g) | 5 — fc is the largest small atom of S} . 

Note that B{g, k) is a subset of B{g) and that B{g, k) is empty whenever k > g/2 by 
lemmajsj In order to describe B{g, k) we are led to the next definition. An ordered 
pair (L, M) of subsets of (0, k) is admissible if it satisfies two conditions: 

(adl) L C M, and 

(ad2) for every x G M there exists y E L with x + y < k and x + y ^ M. 

Let ^(fc) be the set of all admissible pairs of subsets of (0, k), and let Ak = 
denote the cardinality of this set. The power set V{k,g — k) of the set {k,g — k) 
consists of all subsets of (fc, g — k) and has cardinality 2^"^''^^. 

Theorem 5. For integers g and k with g > 2k > the set B{g, k) is in one-to-one 
correspondence with A{k) x 'P{k,g ~ k). In particular, the cardinality of B{g, k) 
equals Ak 29-2fc-i, 
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M+g-k 

/ \ 

-* — \ — \ — \ — \ — \ — * * 
g-k g g+1 



Figure 2. The numerical set S{L, M, P) 



Proof. Suppose that (i, M) £ A{k) and P e ^(fc, g-k). Then 

S{L, M,P)^NgULUPU{g~k}u{g-k + M) (1) 

is a numerical set with Frobenius number g. (See figure [2]) Since g ~ k + L C 
g — k + M by (adl) and each nonzero element of S{L, M, P) — L is larger than k, 
5 — fc is a small atom for S{L, M, P). Suppose that x € {g — k, g) Ci S{L, AI, P). 
Then x — g + k € Ad and by (ad2) there is an integer y G L such that y + x < g 
and y + X ^ g — x + M . This shows that y + x ^ S{L, M, P) and that x is not 
an atom for S{L, M, P). Thus 5 — fc is the largest small atom for S{L, M, P) and 
(L, Ad, P) S{L, M, P) describes a function B from ^(fc) x 7'(fc, g-k) into B{g, fc). 

Now assume that S E B{g,k), and define Ls C (0, fc), A/5 C (0,fc) and P5 C 
{g- k,g) by 

Ls = Sn{0,k), A4s^k-g+{Snig-k,g)), Ps ^ S H (k, g + k) . 

Since g — fc is a small atom of S then £ + g — k E S r\{g — k, g) for each £ E Lg, which 
implies that £ E Ads and that is C A/5. Suppose that x E Ms. Then g — k + x is 
an element of S but not an atom of S (since g — fc is the largest small atom of S) 
and so there exists y E S with g — k + x + y < g and g — k + x + y ^ S. It follows that 
x + y < k, y E Ls and x + y ^ Ads. Thus the pair {Ls, Ads) satisfies (adl) and (ad2) 
and {Ls,Ads) E A{k). Let <I> be the function from B{g,k) to A{k) x 'P{k,g - fc) 
given by 5 1-^ (L5, A/5, P5). The proof is completed by observing that 9 and <i> are 
inverses of each other. □ 

By lemma [s] the set B{g) can be expressed as the disjoint union of the sets 
B{g, fc) where fc ranges from 1 to [(.g — 1)/2J . Thus, the cardinality of B{g) is the 
sum of the cardinaUties of B{g, fc) for 1 < fc < \_{g — \) /2J . 

L(9-1)/2J 

Corollary 6. For each positive integer g, f3g = AkA^^ . 

fc=i 

Proof. Because of theorem [5] and the comment above, we have 

M L(9-1)/2J ,s| L{9-1)/2J . oo-2fe-l L(9-1)/2J 

Pa - 29-1 - 29-1 - 29-1 - -^fe'i • 

fe=l fc=l fe=l 

□ 

By corollary |6] the sequence {/3g} is increasing, and being bounded above by 1, 
it must have a limit 

00 

-fc 



/3oo = hm /3g = V Afc4 



k=l 
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As a consequence the sequence {7^} = {1 — /?g} is decreasing with hmit 

7oo = hm 7g = 1 - /3oo ■ 

By the next lemma, it is also possible to express 700 as the sum of a positive series 
7oo=Er=i(3'=-'-^fc)4-'=. 

Lemma 7. For each integer k > 0, 2L('^~i)/2J < Ak < 3'''^. Moreover 72^-1 — 700 
is positive and 72^-1 - 7oo = ^00 - P2k-i < (3/4)''"\ 

Proof. Let L be an arbitrary nonempty subset of (0, [(fc+l)/2j) C (0,fc) with 
maximal element i. For any element x^L, i+x<2£<k and £ + x ^ L. This 
shows that {L, L) is an admissible pair of subsets of (0, k). Since there are 2L('^^i)/^J 
distinct subsets of (0, [(fc + 1)/2J), this verifies the inequality 2L('=-i)/2J < Ak- 

Suppose (L, M) is an element of A{k). Then for each x G (0, A:) we have three 
distinct possibilities: (1) x ^ M, (2) x e Af and x ^ L, or (3) x ^ L. Therefore 
there are S'^"^ pairs of subsets {L,M) in (0,fc) which satisfy (adl), and it follows 
that Ak < 3*^^^. Now by definition and corollary [g] 

00 00 

72fc-i - 7oo = /?oo - /32fc-i = I] < ^ 4-^ = (3/4)'=-i . 

i— /c i—k 

□ 



Notice that (0, 0) is the only ordered pair of subsets of (0, 1) = 0, and as it 
is admissible, this shows that Ai = 1. Among ordered pairs of subsets of (0,2), 
condition (adl) fails for ({1}, 0) and condition (ad2) fails for (0, {!}) while the two 
remaining ordered pairs (0,0) and ({1},{1}) are in -4(2), and so A2 — 2. With 
lemma |7] and these values of Ai and A2, 

/S„<A + (3/4)=.(i + A) + A.l|, 

which shows that both 700 and /3oo are strictly between and 1 . Using this approach 
with the more extensive data compiled in table [!} we see that = .510538 . . . 
approximates Poo to within (3/4)-'^^ = .0100226.... Taking midpoints gives the 
approximation /3oo ~ .515549 accurate to within .005011, and subtracting from 1 
gives 7oo ~ .484451 with the same degree of accuracy. This approximation can be 
rephrased as saying that \G{g) \ ~ .484451 x 2^~^ for large values of 5. 

If (L, M) is an admissible ordered pair of subsets of (0, k) and L' and M' 
are subsets satisfying L C L' C M' C M then {L',M') is also admissible. In 
particular, both {L,L) and (M, M) are elements of A{k) whenever {L,M) g A{k). 
The computer routine that was used to generate the data in table [T] starts by first 
determining the collection of subsets L C (0, k) for which {L, L) is admissible. 
(The cardinality of this collection is denoted by A'f, in the table. These numbers 
are important in their own right as we shall explain in the next section.) The 
routine then isolates nested pairs of sets in this collection and tests only these pairs 
for condition (ad2). Even with this, the algorithm has exponential complexity and 
slows down quite rapidly. 

From the results of this section we may draw some further conclusions which 
directly address the Anti-Atom Problem for an arbitrary numerical monoid M. 
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n 


A' 


An 




/32n+l + (3/4)" 




1 


1 


1 


.250000 


1.000000 


- 


2 


2 


2 


.375000 


.937500 


.5000 


3 


3 


3 


.421875 


.843750 


.6667 


4 


6 


8 


.453125 


.769531 


.3750 


5 


10 


18 


.470703 


.708008 


.4444 


6 


20 


50 


.482910 


.660889 


.3600 


7 


37 


135 


.491150 


.624634 


.3704 


8 


74 


385 


.497025 


.597137 


.3506 


9 


140 


1065 


.501087 


.576172 


.3615 


10 


280 


3053 


.503999 


.560312 


.3488 


11 


542 


8701 


.506073 


.548308 


.3509 


12 


1084 


25579 


.507598 


.539274 


.3402 


13 


2118 


73693 


.508696 


.532453 


.3471 


14 


4236 


217718 


.509507 


.527325 


.3385 


15 


8337 


635220 


.510090 


.523462 


.3427 


16 


16674 


1888802 


.510538 


.520561 


.3363 



Table 1 . Bounds for /?; 



Theorem 8. Let M 7^ Ng be a numerical monoid with Frobenius number g and 
let g - k be the largest element of M n (0,5). Then |t/(M)| < Afe4-'=2f-i < 
i (3/4)'' X 29-1. 

Proof. If M 7^ Ng is a numerical monoid in S{g) and g — k is the largest element 
in M n (0,(7) then g — k is the largest small atom of every numerical set S with 
A{S) = M. Thus g{M) C B{g,k) and |g^(M)| < \B{g,k)\ = Ak23-^^-^. The last 
inequality follows from lemma [7] □ 

The value of k in theorem [s] satisfies < fc < g/2. Since | (3/4)'' x 2^'^ < 
.25 X 2^-1 is less than .484451 x 2^-^ for all k, we see that among all monoids in 
S{g) the one with largest anti-atom set is always Ng (which is not too surprising 
since G{^g) contains more than 48% of the elements of S{gy). 

As k increases from 1 to [(g — 1)/2J the cardinality of B{g^ k) decreases but the 
number of monoids in B{g^ k) decreases as well. For example, B{g^ 1) contains all of 
the symmetric and pseudosymmetric monoids in S{g), while, at the other extreme, 
^{g^ [{9 ~ 1)/2J) contains only one monoid Dg which is defined by 

Dg = NgU{L(g + 2)/2j}. (2) 



Corollary 9. For each nonnegative integer n, we have \G(D2n+i)\ = An and 

\g{D2n+2)\=2An. 

Proof. It is not difficult to show that Dg is the only monoid in B{g, [{g — 1)/2J). 
(If M is a monoid in B{g, [{g - 1)/2J) then M n ([(5 + 2)/2j ,g) = 0.) Thus Dg is 
the atom monoid of every numerical set in B{g, [{g — 1)/2J), and this implies that 
|g(Dg)| = \Big, [{g - 1)/2J)| = AL(9-i)/2j2^"'L(9-i)/2j-i^ from which the corollary 
follows. □ 
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When k is less than [{g — 1)/2J the set B{g, k) will always contain at least two 
distinct numerical monoids (for example, Ng U {g — k} and NgU {g — k — 1, g — k}). 
Thus Dg is the only monoid in S{g) for which the first inequality of theorem [S] is 
sharp. 

The GENERATING FUNCTION FOR {7^} 

For each integer fc > let A{ky denote the collection of all subsets L C (0, fc) 
for which {L, L) is admissible. Thus A{k)' consists of those subsets L which satisfy 
the condition that for each x G L there is y G L such that x + y < k and x + y ^ L. 
The cardinality of ^(fc)' will be denoted by A'j,. 

Theorem 10. There is a one-to-one correspondence between G{g) and A{gy . In 
particular, \g{g) \ = A'g and jg = A'g /23-\ 

Proof. For each L E A{gy consider the numerical set Ng U L G S{g). If x £ L 
then there is an integer y £ L such that Q < x -\- y < g and x -\- y ^ L, which 
implies that a; + y ^ Ng U L. This shows that x is not an atom of Ng U L, and that 
Ng U -L has no small atoms. Thus the correspondence L 1-^ Ng U L is a function from 
^1(17)' to Q{g), and clearly this function is injective. Now suppose S G G{g) and let 
X £ S n (0,g). Since S has no small atoms, there is an integer y £ S such that 
X -\- y ^ S. Thus X + y < g (since the Frobenius number of S is g), y £ S D {0,g) 
and X -\- y ^ S D {0,g). It follows that S D {0,g) is an element of A{g)', and the 
function L 1— Ng U L is surjective. □ 

Theorem 11. For each fc > 1, A'^^k = 2 A'^^k-i ^""^ ^2fc+i = 2 A'^^. - A^. 

Proof. The first equation follows immediately from lemma |4] and theorem |10| For 
the second equation we have 

k 

A'2k+i = Ie(2fc + 1)1 = \S{2k + 1)1 - \B{2k + 1)1 = 22^- - 5]|e(2fc + 1, t)\, 
and by theorem [5] 

k k 

22fc _ Y}B{2k + 1, ^1 = 4'= - ^ 4!"-^ . 
1=1 1=1 
A similar computation shows that 

1=1 

Combining these gives AJjfe+i ^ 2^42^. = —A^. □ 

The set Q{2k + 1) can be constructed from Q{2k — 1) by a process in which 
each element of Q{2k — 1) will spawn either three or four elements of Q{2k + 1) as 
follows. If 5 e Q{2k — 1) and Q is one of the four subsets of {fc, fc + 1} then let 
S{Q) £ S{2k+\) be given by 

S{Q) - N2fe+i u (5 n [1, fc - 1]) U Q U (2 + (5 n [fc, 2fc - 2])) . (3) 

Since S has no small atoms, S{Q) will not have any small atoms larger than fc + 1. 
Furthermore, if Q is one of 0, {fc} or {fc, fc + 1} then fc + 1 is not an atom of S{Q), 
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M 



* * * 

2k-1 



k-1 k 



Q 



M+2 



S{Q) 



* 

2k+1 



k-1 k+2 



Figure 3. The numerical set S{Q) 



and S{Q) is an element of Q{2k + 1). But if Q equals the singleton set {k + 1} 
then sometimes fc + 1 will be an atom for S{Q), in which case S{Q) is not an 
element of Q{2k + 1). From this we see that \G{2k + 1)| = ^2fc+i ^'^^^ times 
\Q{2k — 1)1 = ^2fc-i nihius the number of elements of Q{2k — 1) which spawn only 
three elements of Q{2k + 1), and, since ^'2^.4.1 = ^A'2k-i ~ by theorem 11 the 
number of elements of G{2k — 1) which spawn only three elements of G{2k + 1) 
equals Afc0 As a result of these comments we can view the union of all the sets 
Q{2k + 1) as the vertices of a downward opening rooted tree in which each vertex 
is directly above the 3 or 4 vertices that it spawns, as pictured in figure |4j In the 
illustration the vertex labeled by a 2 x /c matrix a = ( t?2fc £12" -1 ••• ""+1 ) with entries 
in Z2 corresponds to the numerical set 



S{a) = N2k+i U {i I a, = 1} 

in Q(2k + 1). Although we will not use it here, one can specify conditions on the 
matrix a which guarantee that S{a) is in Q(2k + 1): Call the 2 x fc matrix a 
quadrivalent if there is an integer i with 1 < ^ < fc such that the ith column of a is 
(1) or (1) and the {k+l-£)th column is ( J ) or ( )(] Then S{a) G g{2k+l) if and 
only if whenever the ith column of a equals ( ^ ) then the 2 x (i — 1) submatrix of 
a to the left of that column is quadrivalent. We also note that as one moves down 
the tree the ratio Afe/A2j,_|_]^ of the number of vertices at level 2fc + 1 which spawn 
three vertices by the total number of vertices at that level limits to 0. Indeed Ak is 
bounded above by the number of 2 x fc matrices a which are not quadrivalent and 
that number is easily seen to equal 3'"'. Thus 

Ak Ak Ak ^ 1 /3 ' 



^2fe+i l^(2fc+l)| 72fc+i4'= - 72fc+i V4 

and the latter limits to 0. 

Let g{z) and f{z) be the analytic functions defined by 
00 00 
g{z) = J2 Akz' and f{z) - ^ A'.z'^ . 

k=l k=l 



^This can also be verified by a combinatorial argument. If (L,M) £ A{k) tfien {L',L') is 
an element of A(2k — 1)' = Q{2k — 1) which spawns only three elements of Q(2k + 1), where 
L' = L U (M + A: - 1). 

^The reason for this terminology is that if S{a) is an element of Q(2k + 1), then the matrix a 
is quadrivalent if and only if S{a) spawns four elements of Q{2k + 3). 
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Figure 4. Rooted tree for [j {g{2k + l)\keN}. 



Corollary 12. The functions f{z) and g{z) satisfy the relation 

(2z-l)/(z) = z(g(z2)_l) . 

Proof. Using theorem [TT] we have 

oo oo oo 

k=l fc=l 

Also Ak = ^A'^k-i ^ ^2k+i by theorem [TT] and 

oo 

3(z^) = E (4^2^-1 -^Wi)^'' 

oo ^ oo 

^4zE^k.-i-^'-'-^E^2.-i^ 



fe=l 



fe=l 



^ k=l ^ 



where the last equality follows from Q. 



(4) 



□ 



Corollary 13. The analytic function f{z) has a singularity at z — 1/2, and its 
radius of convergence at the origin equals 1/2. Other than z = 1/2, the singularities 
of f{z) coincide with those of g{z'^) and f{z){z — 1/2) is continuous on \z\ < 1/2. 

Proof. Since J2kLi ^k^^'' sums to (iao — .9(1/4), Y^^=i ^kz'^ converges absolutely 
and f{z) is continuous on \z\ < 1/4. Note that /(z) has no singularities in \z\ < 1/2 
because g{z) has none in \z\ < 1/4. The remainder of the proof follows from 
corollary [12] □ 

By definition the integers A'^, and A/, satisfy < A'j. < A^ for all fc > 0, 



and by corollary 13 we know that the series X^fcLi ^'k^'^ diverges when z — 1/2. 
Therefore J2T=i -^kZ^ must also diverge for z = 1/2 by the comparison test, and 
this shows that the radius of convergence for g{z) at the origin is between 1/4 
and 1/2. To find the precise value, the ratio test would lead one to examine the 
sequence An-i/An- Empirical evidence from the last column of table [T] perhaps 
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suggests that this sequence has a Hmit infimum larger than 1 /4, but we have not 
been able to ascertain this. So the determination of the radius of convergence of 
g{z) at the origin remains open. Notice that if this radius of convergence is larger 
than 1/4 then g{z'^) is analytic in a neighborhood of z = 1/2 and 

^Hm^ {z - l/2)/(z) = ^lim^ z {g{z^) - l) /2 = (/3^ - l)/4 = -7oo/4, 
which would imply that f{z) has a simple pole with residue —700 /4 at z = 1/2. 
Now the generating function for the sequence {'jk} is 

CO 

Hz) = ^ 7fc 
fe=l 

which satisfies 

= E ^ - 2/ (z/2) = (5(^V4) - 1) . 

By corollary [13] the radius of convergence of this series equals 1 . If the radius of 
convergence at the origin for g{z) is larger than 1/4 then (;(z^/4) — 1 is analytic in 
a disk with radius larger than 1 centered at the origin, and h{z) has a simple pole 
at z = 1 whose residue is —700 (and this is the only pole on \z\ = 1). 



Symmetric numerical sets with no small atoms 

A numerical set S with Frobenius number g is negative semisymmetric provided 
that g — X ^ S whenever x G S. For a positive integer g, define 

S'^{g) — {S ^ ^{9) I S is maximal negative semisymmetric in S{g)} 

where maximality is measured with respect to subset inclusion. Then S'' (g) consists 
of the symmetric numerical sets in S{g) when g is odd, and the pseudosymmetric 
numerical sets in S{g) when g is even. Each element of S'^{g) is the union of 
with a subset of {0,g) — {5/2} that is carried onto its complement hy x 1^ g — x, 
and hence the cardinality of S'^{g) equals 2L(9^i)/2J ^ 

The set S'^{g) is partitioned into two subsets 

g'^(g) ^ {S e S^ig) I 5* has no small atoms} = Gig) nS'^ig) 

and 

B^ig) ^ {5* e 5 '"(5) I 5 has at least one small atom} = S(.g) nS^ig) . 
We define 

^,_\B-{g)\ |6-(.9)| 



\S-^ig)\ 2L(9-i)/2J 
and 

|g-(g)| ^ |g-(g)| 
\S''{g)\ 2L(9-i)/2J' 

The next lemma describes a direct connection between symmetric and pseudosym- 
metric numerical sets. 
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Lemma 14. The correspondence S S'q where 

5*0 = (^S* n [0, n - 1]) U (l + n [71, oo)) 

dehnes a bijection from S {2n—l) to5°'(2n) which carries Q {2n—l) ontoQ"'{2n). 
Therefore j^^,^ = and Pi,^ = Pi,^,,. 

Proof. Notice that is the numerical set S'^ where e = as defined in lemma |4] 
The proof follows upon observing that S is symmetric if and only if 5*0 is pseu- 
dosymmetric. □ 

For integers g and k with g > 2fc > 0, let 

B^ig^k) ^ {S e B'^ig) | g - fc is the largest small atom of S} 

(i.e. B'^{g,k) = B{g,k) n S'='{g)). Then B^ig) is the disjoint union of the sets 
B'^{g,k) as k ranges between 1 and [{g — 1)/2J by lemma [s] In order to describe 
B'^{g, k) we are led to the next definitions. 

Let M be a subset of (0, k). Define 

M+ ^ {m e M \ k ~ m e M} , 
M_ = {x e {0,k) \ X M and k - x M} , 

and 

M* ^{x£{0,k)\k~xiM} 
With these definitions, observe that Af* = (M — Af+) U 

A subset M C (0, k) is called a-admissible if it satisfies the two conditions: 
(cr-adl) M- = 0, and 

((T-ad2) for each x £ M there is y £ M* with x + y < k and x + y ^ M. 

Let A'^{k) be the set of all tr-admissible subsets of (0, k) and let denote the 
cardinality of this set. Also, for integers k and g with g > 2k > let V"^ {k, g — k) 
be the collection of all subsets of {k,g — k) — {g/2} that are carried onto their 
complement by the refiection x t—^ g — x. The cardinality of this collection is 
2L(s-2fe-i)/2j ^ jiecall the definition of the numerical set S{L, M, P) as 

S{L, M, P) = Ng U L U P U {g - fc} U {M + g - k) 
from equation ([TJ in theorem [sj 

Theorem 15. Let g and k be integers with g > 2k > 0. The correspondence 
{M,P) S{M*,M,P) defines a bijection from A" {k)xV''{k,g-k) onto B " {g , k) . 
In particular, the cardinality ofB^^g, k) is 2L(9-2fe-i)/2j ^ 

Proof. Let(M,P) e A'' {k)xV (k, g-k). An integer a; e (0, fc) is an element of Af* 
if and only if k — x ^ M, which is equivalent to asserting that g — x ^ M + g — k. 
This together with the fact that P is an element of 'P''{k,g — k) implies that 
S{M*,M,P) e S^ig). Since Af_ = 0, M* = Af - Ar+ C AT. E x £ M* then 
X G M and g — k + x E M + g — k C S{M* , M, P) , and if x is an integer larger than 
k then g — k + x > g and g ~ k + x E S{M* , M, P). This shows that 5 — A: is an 
atom for S{M*,M,P). An element of S{M*,M,P) in the interval {g - k,k) has 
the form x + g—k for some x G M. By the definition of ^"^(fc) there is y E M* with 



This last definition is closely related to the definition of the dual S* of a numerical set S. If 

5 = Nfc U Af 6 S{k) then 5* = Nfe U M*. 
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y < k — x and x + y ^ M. Thus y+(x + g — k) is an element of {g~k,g) which is not 
an element of Af + 5 — fc, and x + g — k is not an atom of S{M* , M, P). It follows 
that 5 - fc is the largest small atom of S{M*,M, P) and S{M*,M, P) eB''{g,k). 

To complete the proof it only remains to show that each numerical set in 
B'^{g,k) equals S{M*,M,P) for some {M,P) £ A"" {k) x V {k,g - k). Let T G 
B''{g,k), and set M = (Tn (g- fc,.g)) - .g + A: C (0,fc) and P ^ T D {k,g-k). If 
X S T n (0, fc) then x + g — k E T D {g — k, g), since g — fc is an atom of T, and 
hence x e Af. Note further that fc — a; is not in T because if it were then both 
k — x and g — {k — x) — x + g — k would be elements of T contradicting the negative 
semisymmetry of T. This shows that x £ M* and Tn {0,k) C A/*. Moreover if 
X £ M* then k — x ^ M which means that g — k+{k — x)=g — x^T and that 
X £ T since T is maximally negative semisymmetric. Thus M* — T (1 (0, k) and 
^(Af*, Af, P) = T. Clearly P £ V''{k,g- k) so to complete the proof it must be 
shown that M £ A''{k). If A/_ ^ then there is a; G (0, k) such that x ^ M and 
k — X ^ M, and then {g — k) + k — x = g — x^T and x ^ T which contradicts the 
maximality of T. This verifies that Af_ is empty. If a; G M then 5 — fc + a; is an 
element oiTr\{0, g) larger than g — k so that g — k + x is not an atom of T since g — x 
is the largest atom of T. It follows that there is an element y £ T with y < k — x 
such that y + g — fc + a: ^ T (note that ?/ cannot equal k — x because otherwise both 
k — x and g — k + x would be elements of T contradicting the assumption that T 
is negative semisymmetric), and M £ A'^{k). □ 

By lemma [3] and the theorem we have 

L(9-1)/2J L(9-1)/2J 
k=l k=l 

and dividing by 2L(9^i)/2J produces the next result. 

L(5-1)/2J 

Corollary 16. For each g > 0, ^ ^ 2-^. □ 

Thus {/Sg*^} is an increasing sequence, and it has a limit 

00 

/3^ = E^^"2-'= . 

It follows immediately that {7/} is a decreasing sequence which converges to = 
1-/3°' 

Moo ■ 

Corollary 17. For each positive integer n, A^^ < 3L("~'^)/^J and 

I2n-1 loo — Poo P2n-1 — 

Proof. Let M C (0,n) be an clement of A'^{n). Suppose that n = 2fc + 1 is odd. 
Then (0, n) is partitioned into k doubletons {i, n — i} where 1 < i < fc. Since 
M- = the intersection of {i,n — i} with M must be nonempty, and so there are 
three possibilities for each of these intersections. Notice also that n—1 cannot be an 
element of M by condition ((T-ad2) , and the intersection of Af with the doubleton 
{l,n — 1} must be {!}. Thus there are at most 3'^^^ = 3L("^'^)/^J possibilities 
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for M. When rt = 2fc is even, (0,ri) can be partitioned into (fc — 1) doubletons 
{i,n — i} where 1 < i < fc — 1 and a singleton {k}. The intersection of M with 
{k} must equal {k} since M_ =0. By a similar argument as before there are at 
most 3'=-2 = 3L("-3)/2j possibilities for M. Now (3^ - I3^^_^ = EfcL„ ^fc" 2"'' < 
^-^3L(.-3)/2j2-.-<(^)"-\ □ 

Some values of l^-^n-i ^i'"^ given in table [2j Note that /3g'^ — .76356 . . . ap- 
proximates f3^ to within (\/3/2)'^^ = .0115731... by corollary 17 Subtracting 
from 1 gives = .23644. . ., which approximates to within .0115731. Taking 
midpoints gives the approximation 7^ « .230653 accurate to within .00579. 



n 


^2n-l 






Rn 




n 




A 




Rn 


1 


1 


1 





1 




17 


16194 


330 


.75290 


.711 


2 


1 





.5 






18 


32058 


206 


.75542 


.744 


3 


2 


1 


.5 


1 




19 


63910 


888 


.75620 


.700 


4 


3 





.625 






20 


126932 


612 


.75790 


.725 


5 


6 


2 


.625 


.871 




21 


253252 


2571 


.75848 


.688 


6 


10 





.6875 


00 




22 


503933 


1810 


.75971 


.711 


7 


20 


3 


.6875 


.855 




23 


1006056 


7274 


.76014 


.679 


8 


37 


1 


.71094 


1 




24 


2004838 


5552 


.76100 


.698 


9 


73 


7 


.71484 


.806 




25 


4004124 


21099 


.76134 


.671 


10 


139 


3 


.72852 


.896 




26 


7987149 


16334 


.76196 


.689 


11 


275 


17 


.73145 


.773 




27 


15957964 


61252 


.76221 


.665 


12 


533 


7 


.73975 


.850 




28 


31854676 


49025 


.76266 


.680 


13 


1059 


43 


.74146 


.749 




29 


63660327 


179239 


.76285 


.659 


14 


2075 


24 


.74670 


.797 




30 


127141415 


146048 


.76318 


.673 


15 


4126 


118 


.74817 


.728 




31 


254136782 


523455 


.76332 


.654 


16 


8134 


74 


.75177 


.764 




32 


507750109 


440980 


.76356 


.666 



Table 2. Approximating 13^, where i?„ = 1/ 



For a numerical monoid M G S{g) let 

G^iM) = g{M) nS^ig) = {Se S^ig) \ A{S) = M} . 

Note that M will not be an element of Q°'{M) unless M is symmetric or pseu- 
dosymmetric, and that Q'^(M) may be empty. If M g S{g) is a numerical monoid 
in B{g,k) (which means that g — k is the largest integer in M n (0,17)) then 
G^iM) (ZB''{g,k). Therefore 

\g%M)\ < fc)| = 2L(f-2fe-i)/2j <^(^\ X 2L(s-i)/2J 

3v3 \ 2 y 

by corollary 1 17| and this is the symmetric analogue of the inequality in theoremjS] In 
particular we see that |tJ'^(Ng)| = |5'^(.g)| is larger than |5°'(M)| for every numerical 
monoid M G S{g) other than Ng. Taking k = [{g — 1)/2J in the above shows that 
\g''{B)2k+i)\ = \B''{g,k)\ = and that IG" iB)2k+2)\ = ^A^, where D„ is defined 
in equation ([2]). 
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The GENERATING FUNCTION FOR {7^^"} 

We define A'^{ky to be the subset of A'^{k) consisting of all cr-adniissible sets 
M C (0, k) for which M+ has at most one element. If M G A'^{ky then there are 
two possibilities: either k is odd and = (because Af+ has an even number of 
elements whenever k is odd) or k is even and = {k/2} (because fc/2 must be 
in M_ or M_|_ whenever k is even, and M_ = 0). The cardinality oi A'^{k)' will be 
denoted hy A^'. 

Theorem 18. There is a one-to-one correspondence between the sets G'^{g) and 
A^ig)'. In particular, IG^ig)] = A^' and 7/ = yl^^' 2-L(ff-i)/2J . 

Proof. Given M e A^ig)' define 5 = Ng U M* . If g is odd then M_ = M+ = 
so that M* = M and S is symmetric. If g is even then M_ = 0, M_|_ — {g/2} 
and M* — M — {5/2} which implies that S is pseudosymmetric. In either case, S 
is an element of S'^{g) and we can write M* = M — {g/2}. Suppose s is a large 
small atom in S. Then s € M* C M and by (cr-ad2) there exists y G M* such that 
s -\- y < g and s -\- y ^ M. Note that s -\- y ^ M — {17/2} = M* since s is large. 
Thus s + y ^ 5 in contradiction of the assumption that s £ ^(5*). This shows that 
M 1-^ is an injective function from A'^{g)' to Q'^{g). If S G G'^{g) then it is not 
hard to check that S = NgUM* where M = {Sr\{0,g))* e A^ig)', completing the 
proof. □ 

Theorem 19. For each fc > 1, A^,! = ^2k-i and A^,!^^ = 2A:^^ - A^ . 

Proof. The first statement follows from lemma [14] and theorem [18] For the second 
we have 

k 

A2k+i - |e"(2A: + 1)1 = |5"(2A: + 1)1 - |6"(2A: + 1)1 = 2'= - ^ A/2'^-^ 

1=1 

using equation and A21! — 2*^^^ — X^fci A^2^~^~^ by a similar computation. 
Combining these two equations gives 2A^j.' — j4^j,'_j^ = A^ . □ 

As in figure[4]we may view IJ {5 {2k + 1 ) | A: e N} as the vertices of a downward 
opening rooted tree. Here each element oiQ" {2k~'\) (represented by a vertex in the 
fcth level of the tree) will spawn either one or two elements of Q'^{2k-\-l) (represented 
by vertices at the (fc + l)st level). If S" S Q'^{2k - 1) and Q is either {fc} or {fc + 1} 
then let S{Q) be defined by equation ([S]) and observe that S{Q) £ S^{2k + 1). As 
before S{Q) has no small atoms when Q = {fc} but fc + 1 may be a small atom for 
S{Q) when Q = {fc + 1}. From this we see that |^?'^(2fc + 1)| = ^sV+i is two times 
\Q'^{2k — 1)1 = A2-f!_i minus the number of elements of G'^{2k — 1) which spawn 
only one element of Q"{2k +1), and it follows from theorem [19] that the number 
of elements of Q"'{2k — 1) which spawn only one element of Q'^{2k + 1) equals A^. 
Figure jS] shows the first few levels of the rooted tree. In this illustration a labeling 
sequence a = (ai, . . . , ak) with entries in Z2 represents the same numerical set in 
^ °'(2fc+ l)as the 2 x fc matrix ( ... ) represented in figure|4j where a* = 1 — a^. 
Call the sequence a bivalent if there is an integer £ with 1 < ^ < fc such that 
ai = ak+i-i = 1. Then a represents an element of Q'^(2k + 1) if and only if 
(ai, . . . , ai-i) is bivalent whenever ai = 0. 
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Figure 5. Rooted tree for [J {g''{2k + 1) \ k e N}. 



If we associate a = (ai, . . . , a^) with the finite set F{a) — {0} U | = 1} 
then a is bivalent if and only if k + 1 G F{a) + F{a). Furthermore, a represents 
an element of Q'^{2k + 1) precisely when F{a) is an 'additive 2-basis for k' (which 
means that [0,A:) C F{a) + F{a)). Thus A^flj^^ equals the number of subsets 



of [0, k) which are additive 2-bases for k. By theorem 18 A^^j^i = 1^"^^ which 
is asymptotically equal to "i^2^ where 7^ w .230653 as described above. The 
study of finite additive 2-bases for k has a long history, especially in relation to the 
determination of bounds for the smallest cardinality of such bases. The introduction 
of |GNj has a nice overview of this. The first 19 terms of the sequence 
have been posted at [S] by M. Torelli as sequence number A008929. The paper [T] 
describes some related sequences. (In that paper a finite additive 2-basis is called 
a (finite) 'Goldbach sequence'.) 

Let g'^{z) and f^iz) be the analytic functions defined by 

00 oo 

5-(z) = ^ A.-z'^ and r(z) = ^ 
fe=i fe=i 

Corollary 20. The functions f'^{z) and g'^{z) satisfy the relation 
{2z'-l)r{z)^z{z+l){g'^{z')~l) . 



Proof. First observe that 



k=l 



fe=l 



z + i)J2a 

fe=l 



_2fe-l 



(6) 
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by theorem 19 From the same theorem, Aj^ = 2A2f!_i — ^2fe+i ^^"^ 



oo 



z 

k=l k=l 



□ 



Corollary 21. The analytic function f'^{z) has singularities at z = ±l/\/2 and 
its radius of convergence at the origin equals 1/V2. Except for z — ±l/\/2 and 
possibly for z — —1, the singularities of f"'{z) coincide with those of g'^{z'^) and 
/°'(z)(2z^ — 1) is continuous on the closed disk \z\ < 1/V2. 

Proof. Since g'^(l/2) = f3^ the power series X^fcLi ^k^'^ converges absolutely and 



g'^{z) is continuous on \z\ < 1/2. By corollary 20 /"^{z) has singularities at z = 
±1/a/2 and /"{z) has radius of convergence l/v2at the origin. The last property 
also follows immediately from corollary [20] □ 

Since < A^' < , the series J^'kLi (VV2) diverges by comparison with 
'^'kLi ^k' (l/v^)'^, and so the radius of convergence of g"{z) at the origin must 
be between 1 /2 and 1 /\/2- The root test would equate this radius of convergence 
with the limit infimum of i?„ = 1/ y^A^. This value seems to be larger than 1/2 
by the data in the last column of table [2| but we have not been able to ascertain 
this. 

Let h'^{z) = X^fcLi 7fc •^'^ t)e the generating function for {7^}. Using equation ^ 
and corollary [20| we have 

°o ACT' ^ Aal °° A<yl 

ucr/^ _ ^k ^k ^2k-l 2k-l , ^2fc 2k 

k=l k=l k=l 

CO / \ 2/c— 1 

= ^/2(z + l) E ^2^-1 j = 2 (z + 1/z + ^/2) T (V^/2) 
(.9^ - 1) . 

Therefore h"'{z) has radius of convergence 1 at the origin and has a singularity 
at z = 1. If the radius of convergence of g'^(z) at the origin is larger than 1/2 
then z = 1 is the only singularity of h"(z) inside a circle with radius larger than 1 
centered at the origin and h^^z) would have a simple pole at z = 1 with residue 

lim(z - l)h^z) = hm z(g'^(zV2) - 1) - g-(l/2) - 1 = -7^ . 
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